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EVOLUTION  OF  A  VOLTACE-TIME  MOtEL  OF  THERMAL  BATTERIES 

by 

J.  Knight 


SUMMARY 


A  temperature-time  model  of  thermal  batteriea  hae  almost  been  completed  and 
validated  age.inat  experimental  data.  This  Memorandum  first  summarises  early 
atLempts  to  integrate  a  voltage-time  model  into  thii,  taking  advantage  of  the 
instantaneous  prediction!)  of  temperature)  thermodynamic  potentials,  and  internal 
reeietance  which  the  thermal  model  provides.  It  than  daecribce  how  recent 
refinements  of  the  voltage-time  model  have  led  to  improved  simulation  of  the  dis¬ 
charges  of  a  wide  range  of  sizes  and  types  of  thermal  battery  under  an  equally 
wide  range  of  test  conditions. 

The  semi-empirical  approach  adopted  hai  been  to  provide  a  un-varsally- 
applicable  framework  based  on  logical  concepts  to  cover  various  effects  such  as 
polarisation  etc,  but  with  adjustable  numerical  parameters.  It  is  shown  that  a 
moderately  good  simulation  may  be  obtained  for  the  majority  of  available  dis¬ 
charge  curves,  using  this  one  set  of  equations  and  without  filtering  parameter 
values.  Further  improvements  can  be  obtained  when  parameter  values  are  optimised 
for  one  particular  type  of  battery, 

Comparisons  of  model  simulations  against  a  body  of  experimental  data  have 
pinpointed  remaining  discrepancies  which  will  guide  further  rafinemant  efforts. 


2 


LIST  OF  CONTENTS 


Page 


1  INTRODUCTION  3 

2  DESIGN  AND  PERFORMANCE  DATA  BANK  4 

3  THEORETICAL  CONTENT  OF  THE  VOLTAGE-TIME  MODEL  5 

4  PRELIMINARY  VOLTAGE-TIME  MODELS  6 

5  MARK  II  VOLTAGE-TIME  MODEL  7 

6  MARK  III  VOLTAGE-TIME  MODEL  10 

6.1  Capacity  degradation  II 

6.2  Allowance  for  time-dependant  polariiation  II 

6.3  Improvement!  to  model  procedures,  efficiency,  and  output  14 

6.4  Effect  of  discharge  state  on  internal  resistance  15 

6.5  Effect  of  manufacturing  tolerances  17 

7  TEMPERATURE-LIMITED  BATTERY  DISCHARGES  18 

7.1  Salt  freezing  temperature  in  discharged  batteries  18 

7.2  Ocher  heat  aourcea  19 

8  REMAINING  DISCREPANCIES  21 

9  CONCLUSIONS  22 

Acknowledgment*  23 

Appendix  Estimated  magnitude  of  heat  generation  from  aide  reactiona  25 

Table  I  Details  of  batteries  and  teat  condition*  27 

Ra fa ranee a  28 

Illustration*  Figures  1-15 


Report  documentation  page 


inside  back  cover 


!'  •  - 


{ 

I 

f- 


TM  MS  1163 


0096948 


CONDITIONS  OF  RELEASE 


300B73 


MR  PAUL  A  ROBEY 
OTIC 

AttfllOTIC-FOAC 
Camaron  Slillon-Bldg  6 
Alaaandrla 
VA  81304 1148 
USA 


COPYfllQHT  to) 

1888 

CONTROLLER 
HMSO  LONDON 


DR1CU 


. DR,0  Y 

Raporta  quoted  irt  not  naoaaaarlly  aval  labia  to  memberi  of  tht  public  or  lo  oommarolal 
organ  iMtiori*. 


3 


I  INTRODUCTION 

Thermal  bacterial  were  developed  after  World  War  II  as  the  preferred  power 
source  for  various  functions  on  guided  missiles.  They  are  active  at  temperatures 
of  typically  300°C  and  rely  on  molten  salt  eutectics  such  as  lithium  chloride- 
potassium  chloride  in  order  to  supply  significant  power*.  Thus,  compared  with 
most  other  types  of  bateary,  their  performance  is  more  then  usually  influenced 
by  temperature.  This  might  eevtrely  curtail  their  active  lives  if  the  self 
eutectic  solidifies  on  cooling  before  their  electroactive  materials  are 
exhausted.  Equally,  too  high  a  peak  temperature  during  battery  activation  may 
result  in  premature  battery  failure.  Another  factor  which  is  more  extreme  in 
thermal  battery  discharges  than  in  other  types  is  the  wide  range  of  temperature 
experienced  by  their  cells  throughout  their  active  life.  This  varies  according 
to  application,  ambiant  temperature ,  and  siss  of  battery,  but  might  typically  be 
from  5S0  down  to  A00°C  or  below.  Lastly,  the  heat  which  eventually  warm*  up  the 
battery  outer  container  may  sometimes  pose  a  problem  for  the  immediate  surround¬ 
ing!  of  the  battery. 

From  the  above  considerations,  it  is  teen  that  a  temperature-time  model  of 

thermal  batteries  can  of  itself  be  a  conaidarabla  aid  to  the  thermal  battery 

designer.  Predicting  the  temperature  of  older  battery  designs  using  the  calcium/ 

calcium  chromate  electrochemical  couple  was  difficult  because  of  the  extensive 

chomicsl  aide  reactions  which  occur.  However,  thair  increasing  replacement  by 

lithium/iron  disulphide  design*  reduced  these,  simplifying  the  prediction  of  battury 

temperature  and  performance.  Several  organisation*  have  published  their  efforts 

2 

on  thermal  modelling  of  thermal  batteries,  including  Sandia  Laboratories  ,  Their 
modal  which  predicted  temperature  diatributions  during  the  cooling  phase  was 

extended  when  RAE  simulated  the  large  temperature  gradients  built  up  during  acti- 

•1 

vation  by  a  pyrotechnic  charge  ,  Predictions  were  then  available  of  the  differ¬ 
ing  heating  rates  of  the  various  cell  components,  their  peak  temperatures,  and  of 
the  overall  activation  time  for  different  designs.  More  recently  RAE  effort  has 
focussed  on  validating  the  thermal  model  during  the  subsequent  cooling  phase 
against  experimental  thermocouple  data  generated  from  both  electrically-inactive 
and  -active  batteries.  Quite  good  simulation*  of  battery  temperature  are  now 
available  uaing  this  for  a  wide  range  of  battery  design*  and  electrical  load 
profiles.  This  is  now  regarded  as  almost  complete,  the  major  deficiency  being 
that  further  study  would  be  desirable  of  the  extent  of  chemical  side-reactions 
such  as  cathode  pyrolysis  within  the  thermal  battery  environment,  as  opposed  to 
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open  systems  with  inert-gas  flow,  such  as  the  pan  of  a  thermogravimetric 
analyser"*.  Further  discussion  of  exothermic  side-reactions,  as  they  affect  the 
voltage-time  model,  is  given  in  sections  6.1  and  7.2. 

The  final  aim  of  thermal  battery  modelling  studies  at  RAE  continues  to  be 
the  provision  of  an  integrated  temperature  and  voltage  model  which  is  capable 
of  reasonably  accurate  prediction  for  a  wide  range  of  batteries  and  test  con¬ 
ditions.  This  Memorandum  describes  a  further  stage  in  the  improvement  of  the 
voltage-time  model,  as  well  as  documenting  how  this  stage  has  evolved,  and  the 
approach  adopted.  The  advantage  of  the  model  when  finalised  should  be  that  the 
process  of  optimising  new  designs  for  new  applications  should  be  speeded  up,  with 
an  attendant  cost  saving;  it  should  also  aid  research  on  new  configurations  and 
should  increase  the  confidence  with  which  designs  outside  the  presently-known 
range  might  be  configured.  All  this  can  be  assessed  'on  paper'  before  prototypes 
are  manufactured  and  tested. 

2  DESIGN  AND  PERFORMANCE  DATA  BANK 

In  essence,  the  RAE  model  is  fed  detailed  information  on  the  battery  design, 
and  on  the  chosen  electrical  load  profile  and  other  test  conditions;  it  then 
predicts  battery  performance.  A  body  of  combined  design  and  performance  data  is 
thus  essential  against  which  to  test  the  effect  of  model  changes.  As  the  aim  is 
to  eventually  model  a  wide  range  of  battery  types  and  test  conditions,  this 
data  bank  should  also  cover  the  full  range.  RAE  do  not  manufacture  thermal 
batteries,  not  even  experimental  prototypes,  so  the  experimental  data  quoted  in 
this  Memorandum  is  necessarily  limited  and  has  either  been  reproduced  by  permis¬ 
sion  from  thermal  battery  manufacturers'  own  tests,  or  has  been  obtained  by  RAE 
discharge  tests  on  their  products. 

Nevertheless,  the  11  battery  designs  studied  span  a  range  of  outside  dia¬ 
meters  from  18  to  91  mm,  employ  lithium-iron,  lithium-aluminium  or  lithium-silicon 
as  anode  materials,  but  all  use  iron  disulphide  as  the  cathode.  Likewise,  the 
16  discharge  test  conditions  range  from  ambient  temperatures  of  -45  to  +74°C,  and 
electrical  load  profiles  ranging  from  continuous  constant-current  or  constant- 
resistive  loads  to  periodic  sequences  of  short  pulses.  The  designs  are  not 
necessarily  regular  production  designs,  nor  are  the  load  profiles  necessarily 
representative  of  any  particular  application.  Table  1  gives  some  information  on 
these  battery  tests,  and  indicates  whether  temperature  or  capacity  exhaustion 
were  predicted  to  be  the  life-limiting  factor  under  the  conditions  employed. 
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3  THEORETICAL  CONTENT  OF  THE  VOLTAGE-TIME  MODEL 

Ideally,  computer  models  should  be  developed  entirely  according  to  established 
theories,  so  that  they  might  eventually  have  good  predictive  power  beyond  the 
presently-known  limits,  and  might  considerably  increase  understanding  of  how 
observed  discharges  are  derived.  It  was  a  practical  proposition  to  base  the 
thermal  model  almost  entirely  on  the  well-known  theories  of  heat  transfer,  ther¬ 
modynamics  of  electrochemical  cells  etc,  so  that  an  ab  initio  model  could  be  devel¬ 
oped  in  a  reasonable  time. 

This  would  also  be  the  ideal  for  the  voltage-time  model.  Lawrence 
Berkeley  Laboratory,  California  have  described  one  such  model  of  flooded- 
electrolyte  lithium  alloy-iron  disulphide  cells  operated  isothermally6 .  This  is 
developed  from  their  earlier  model  of  iron  monosulphide  cells7,  and  considers 
dek  ils  of  the  volume  make-up  and  porosity  of  the  electrodes,  diffusion  rates 
within  them,  and  che  effects  of  phase  precipitation  on  performance.  The  impetus 
for  their  work  was  the  optimisation  of  the  large  secondary  well-insulated 
batteries  intended  for  vehicle  propulsion,  rather  than  the  comparatix'ely  short¬ 
lived  smaller  primary  thermal  batteries  which  would  usually  experience  a  more 
variable  temperature  during  their  active  lives. 

g 

At  the  other  extreme,  Shepherd  and  others  have  shown  how  good  fits  to 
observed  discharges  of  several  types  of  battery  may  be  obtained  by  purely  empiri¬ 
cal  equations  with  fitting  parameters  readjusted  for  each  type,  and  without  any 
reference  to  the  thermodynamics  of  the  cell  reaction.  In  deciding  which  approach 
might  be  most  practical  to  adopt  for  thermal  batteries,  the  following  features  of 
thermal  battery  discharges  (some  of  which  are  experimental  observations  and  some 
likely  hypotheses)  might  usefully  be  listed: 

(i)  the  discharge  of  iron  disulphide  cathode  can  be  highly  inhomogeneous^; 

(ii)  experimental  discharge  curves  can  be  significantly  modified  by  addi¬ 
tives  ,  or  even  by  different  particle  sizes  or  batches  of  active 
material ; 

(iii)  the  possibility  exists  of  hidden  leakage  currents  (shorting)  during 
the  active  life  of  thermal  batteries; 

(iv)  the  electrodes  and  separator  are  porous  multi-phase  bodies,  and  dis¬ 
charge,  which  depends  on  migration  of  species  through  them,  could 
well  be  influenced  by  factors  like  aspect  ratio,  void  content6, 
wettability  etc. 
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Although  this  list  is  nut  meant  to  be  exhaustive ,  it  gives  sufficient  indication 
that  an  attempt  to  provide  a  truly  theoretical  model  of  thermal  battery  discharges 
might  not  he  coat-effective,  since  at  least  some  of  the  above  affects  could  be 
ill-dsflnad  and/or  difficult  to  determine. 

It  wee  therefore  considered  beat  to  provide  a  theoretical  framework  for  the 
gross  effects  such  as  Faraday's  Law,  cell  thermodynamics,  end  internal  resistance 
which  are  relatively  simple  to  treat  theoretically,  whilst  complementing  this  with 
a  semi-empirical  treatment  of  effects  such  a*  polarisation.  Where  possible  this 
semi-empirical  portion  is  bated  on  physically  reasonable  ideas  end  the  aim  has 
been  to  eat  up  equations  having  general  applicability  for  all  conditions  met. 

These  equations  contain  numerical  parameters  which  are  set  by  fitting  procedures. 
Ideally,  the  same  numerical  values  would  than  be  valid  ovar  all  conditions,  but 
otherwise  there  should  at  leaat  be  some  category  of  conditions  (for  example 
batteries  from  one  particular  manufacturer,  or  one  particular  combination  of 
pellet  material  compositions)  over  which  reliable  predictions  could  be  made. 

The  extent  to  which  these  aims  have  been  achieved  in  the  latest  model  version  is 
apparent  in  section  6. 

4  PRELIMINARY  VOLTAGE -TIME  MODELS 

The  thermal  and  electrical  behaviour  of  thermal  batteries  are  interdependent 
in  a  number  of  important  respects  i 

(i)  the  passage  of  current  through  the  battery  generates  heat; 

(ii)  the  equilibrium  potentials  of  both  electrodes  are  teroperaturr- 
dependenci 

(ill)  the  active  lifetime  of  the  battery  may  be  linked  to  the  freeeiug 
of  molton  salts; 

(iv)  the  internal  resistance  of  the  battery  depsnds  on  temperature. 

The  aim  of  the  first  voltage  model  was  simply  to  estimate  the  rate  of  heat 
generation  caused  by  current  flow,  as  in  (i)  above,  and  so  to  improve  the  thermal 
model.  The  magnitude  and  offset  of  thia  can  be  very  considerable  for  large 
batteries  discharged  at  high  rates,  even  to  the  extant  of  causing  n  net  tempera¬ 
ture  rise  within  the  cells  as  tha  discharge  proceeds,  as  this  heating  effect  can 
more  then  offset  tha  natural  cooling  of  the  battery  according  to  Newton's  Law. 

In  thie  version  the  internal  resistance  value  was  e  conetant  supplied  by  the 
user. 

In  the  first  major  attempt  to  develop  a  voltage-time  model  per  se,  dis¬ 
charge  curves  from  isothermal  single-cell  tests  were  analysed  and  a  number  of 
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concepts  affecting  the  cell  capacity  anil  voltage  were  inferred  and  expressed  in 
mathematical  form.  These  included  the  variation  of  internal  resistance  with  both 
temperature  and  seate-of-disehsrge,  and  Che  degradation/under-utilisation  of 
theoretical  capacity  as  a  function  of  current  density  above  a  certain  threshold. 
Thermodynamic  potentials  were  calculated  only  for  the  complete  cell,  not  for 
individual  electrodes,  so  it  was  assumed  that  there  wee  always  an  exact  balance 
of  theoretical  capacity  between  the  two.  Whilst  good  simulations  were  initially 
obtained  over  a  narrow  range  of  battery  types,  it  was  soon  clear  that  this  model 
would  not  find  more  ganerel  applicability,  and  was  therefore  discontinued. 

5  MARK  II  VOtTASB^TIHE  MODEL 

The  main  improvement  here  ovor  the  previous  attempts  was  to  base  the  calcu¬ 
lation  on  a  sounder  theoretical  footing  by  treating  each  electrode  individually 

9  j  o 

end  incorporating  literature  aquations  ’  ,  defining  tho  various  electrode 

potentials  as  a  function  of  both  state-of-diaeharge  and  temperature.  ThiB 
automatically  catered  for  battery  designs  in  which  the  electrodes  were  unbalanced 
in  theoretical  capacity.  Apart  from  this,  this  version  essentially  computed  the 
discharge  curve  one  might  expect  from  a  battery  in  which  the  active  ingredients 
were  100%  utilised,  without  any  kinetic  limitation,  and  assuming  no  dissipation 
of  theoretical  capacity  other  then  usefully  through  the  external  load.  As  a 
minor  variant  incorporated  later,  an  ability  to  simulate  short-circuit  paths  was 
provided,  if  a  constant  resistance  value  was  entered.  The  Mark  II  model  was 
therefore  too  optimistic  for  manufactured  batteries,  as  expected,  except  in  the 
important  caee  where  the  test  conditions  imposed  cell  cooling  ss  the  factor 
limiting  the  battery  life.  This  case  was  already  simulated  by  the  existing  ther¬ 
mal  model  and  examples  of  its  effects  have  already  been  given4  end  are  also  shown 
later,  Whilet  being  too  optimistic  for  discharges  which  are  not  limited  by 
temperature,  this  version  was  of  interest  in  comparing  experimental  discharges 
with  the  theoretical  limit,  in  order  to  assess  energy  efficiency  more  clearly. 

The  following  steps  summarise  how  tho  prediction  of  on-load  voltage  is 
derived i 

(i)  the  user  defines  the  dimensions  and  densities  of  cell  component 
pellets  in  his  input  data| 

(ii)  the  user  also  specifies  a  number  of  either  constant-current  or 
constant-resistive  loads  for  a  combination  of  both)  together  with  the 
time  period!  during  which  each  operates; 

(iii)  from  (i)  the  theoretical  capacities  of  both  electrodes  are  computed 
(for  pellets  with  the  minimum  specified  dimensions). 
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(iv)  the  number  of  celle  is  derived  by  counting  Che  number  of  oleetrolyte/ 
cathode  interfaces  in  the  user's  design  grid: 

(v)  adjustments  are  made  to  all  affected  variables  (number  of  cells, 
current,  internal  resistance)  according  to  whether  the  user  chose  a 
quarter-  or  half-section  grid,  or  whether  h#  specified  one  series 
string  of  cells,  or  two  equal  stacks  connected  in  parallel; 

(vi)  the  iterative  process  which  recalculates  each  variable  at  increasing 
times  begins; 

(vii)  the  currant  flowing  through  the  cells  is  computed,  in  cases  (eg 
short-circuiting),  where  it  is  not  identical  to  that  in  the  external 
circuit,  and  for  either  type  of  external  load; 

(viii)  from  (vii),  a  running  tally  is  kept  of  discharged  capacity,  and 
resistive  heating  and  antropic  cooling  are  estimated) 

(ix)  temperature  distributions  are  calculated,  corrected  for  phase 
changes,  and  a  single-value  'cell  temperature'  is  defined  as  the 
mean  of  the  maximum  and  minimum  cathode  temperatures; 

(x)  a  atate-of-discharge  ia  defined  as  the  ratio  of  the  discharged  capa¬ 
city  to  that  theoretically  available  from  a  l-elaet.ron  transfer  in  the 

Q 

computed  amount  of  active  material  in  each  electrode  j 

(xi)  a  thermodynamic  equilibrium  potential  for  each  eleerrode  is  computed 

9  1 0 

from  the  literature  equations  mentioned  above  '  ; 

(xii)  the  cell  and  battery  emfr  are  defined  from  the  algebraic  sum  of 
electrode  potentials) 

(xiii)  numerical  values  are  assigned  to  the  bulk  resistivities  of  all  coll 
component  pellet*.  These  are  a  mixture  of  inferred  values  from  a 

9 

report  together  with  arbitrary  interim  values  pending  better 
measurements.  The  values  used,  in  ohm  cm,  were) 

Pellet  material  Resistivity  Reference 


lithium/iron 

0.3 

lithium/aluminium  or  lithium/silicon 

0.5 

iron  disulphide  cathode  pellets 

1.4 

9 

all  electrolyte  pellets 

2.68 

9 

burnt  pyrotechnic  pellets 

1.0 

metallic  cups  or  discs 

0.1 

(The  variation  of  cathode  resistance  with  stete-of-discharge  quoted  in  a  report9 
was  also  ineorporatsd,  but  was  withdrawn  after  it  resulted  in  poorer  results  for 
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on*  or  two  batteries  simulated .  A  possible  reason  for  this  has  now  been  given, 
as  discussed  in  section  6.4,) 

(xiv)  the  above  electrolyte  resistivity  is  then  modified  if,  after 
inspection  of  the  temperature  of  every  electrolyte  grid  element, 
it  is  found  that  salt  eutectic  would  either  be  within  or  below  the 
freeeing  range  in  the  element,  These  extra  contributions  are  then 
summed  across  each  cell  pellet,  and  then  down  each  series  string  of 
cells,  in  order  to  estimate  the  total  effect  of  salt  freezing,  if 
any,  throughout  the  battery! 

(xv)  an  interim  internal  resistance  for  the  whole  battery  (Irrespective 
of  type  of  design  grid  or  electrical  connection)  is  then  computed 
from  the  resistivity  data  ((xiii),  (xiv))  and  from  the  pellet 
dimensions  the  user  specified  (l).  This  ia  then  corrueted  for 
temperature; 

(xvi)  an  on-load  voltage  is  then  estimated  at  aach  sample  time  from  the 
aimple  formula 

V  -  V'  -  1R  ,  (I) 

where  V  is  the  battery  on-lood  voltage, 

V'  is  the  battery  thermodynamic  equilibrium  potential  (xii), 
I  ie  the  current  flowing  through  the  cells  (vii) 

and  R  is  the  corrected  internal  resistance  (xv)  at  the  cull 
temperature  (ix)i 

(xvii)  to  balance  economy  of  computing  time  with  the  ability  to  calculate 
and  display  pulse  shapes,  a  programme  of  calculation  times  is  set 
up,  in  which  the  interval  between  successive  calculations  starts 
increasing  by  'it  each  time  once  the  large  temperature  gradients 
formed  during  activation  have  subsided.  At  the  approach  to  each 
electrical  load  change  specified  by  the  user,  this  programme  is 
temporarily  interrupted  so  that  a  calculation  may  ba  done  immedi¬ 
ately  on  either  side  of  the  load  change.  Afterwards  the  programme 
of  expanding  calculation  intervals  is  mads  to  resume  from  the  value 
it  reached  just  before  the  load  change  occurred; 

(xviil)  the  on-load  voltage  (as  well  as  temperature  distributions  and  other 
heat  loss  parameters)  are  sampled  at  all  tines  specified  by  the 
user,  and  automatically  on  either  side  of  each  electrical  load 
change,  whether  they  ware  specifically  requeuted  or  not. 
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Thi  effect*  of  applying  the  Mark  XX  model  to  the  discharges  of  seven  dif¬ 
ferent  battery  designs  (Table  1,  Batteries  t-7)  discharged  on  either  continuous  or 
pulsed  load*  are  shorn  in  Figs  I  to  3.  (Except  where  otherwise  stated,  all  solid 
traces  in  the  Figures  represent  experimental  data)  broken  traces  are  model  simu¬ 
lations.)  The  following  features  of  the  simulations  may  be  listed; 

(i)  too  low  a  battery  voltage  is  predicted  immediately  after  activation 
for  eoffla  batteries  (Fig  2,  battery  2;  Fig  3).  This  is  probably  dua 
to  the  well-known  'voltage  spike'  effect  which  the  model  does  not 
allow  fort 

(ii)  after  this  has  gone,  the  model  estimates  the  battery  voltORQ  quite 
wall  early  in  discharge  life  (Figs  2  and  3),  and  where  it  is  pos¬ 
sible  to  check  on  the  magnitude  of  internal  resistance  (Fig  2, 
battery  2;  Fig  3)  the  resiativity  values  quoted  above  provide  quite 
a  close  estimate,  This  is  especially  apparent  for  battery  3  in 
Fig  3.  Thus  the  calculations  of  both  thermodynamic  potentials  and 
internal  resistance  are  reasonably  acmrate  early  in  life; 

(iil)  there  is  a  consistent  tendency  for  the  simulated  battery  voltage  to 
remain  steady  in  later  life,  whereas  in  reality  it  gradually  falls 
away  sooner  than  was  predicted.  The  only  exception  seen  is  one  of 
Che  predicted  discharges  of  battery  6  in  Fig  2  (which  differ  only 
In  teat  temperature),  whichwas  predicted  to  be  tempernturo- 
1  limited  due  to  the  low  ambient  temperature) 

(iv)  for  batteries  subjected  to  regular  but  occasional  pulses  on  top  of 
an  otherwise  low  or  zero  background  rurrenc  (FiR  3),  and  which  are 
predicted  to  be  temperaturu-lim-ted,  the  Mark  IX  modal  provides  for 
an  increased  internal  resistance  as  eventually  observed,  but  in  the 
simulation  this  appears  too  early  in  life  and  increases  too  suddenly 
in  magnitude. 

Similar  trends  were  observed  fur  the  other  few  batteries  in  Table  ).  Fig  1 
shows  how  this  model  naturally  leads  to  the  formation  of  'square'  pulses  rather 
than  the  ones  observed  experimentally.  It  is  also  seen  that  in  moBt  canes,  too 
optimistic  a  duration  would  have  been  predicted  to  any  particular  voltage  cut-off. 
The  next  section  describes  how  some  at  least  of  theae  inadequacies  were  addressed. 

6  MARK  XII  V0LTACE-T1ME  MODEL 

This  enhancement  over  the  Mark  II  version  provides  a  preliminary  treatment 
of  polarisation,  reintroduces  a  simpler  concept  of  capacity  degradation,  makes 
allowance  for  manufacturing  tolerances,  and  revises  the  presentation  of  results, 
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6.1  Capacity  degradation 

Comparison  between  the  theoretical  capacity  of  batteries,  their  observed 
duration  to  exhaustion,  and  the  amount  of  residual  lithium  metal  estimated  from 

post-mortem  analysis  after  the  battery  haa  cooled  shove  that  the  active  materials 
can  be  partially  consumed  at  the  same  time  at  they  provide  useful  electrical  out¬ 
put'  1 ,  There  ia  also  evidence  from  a  very  limited  number  of  tests  that  the  degra¬ 
dation  rata  can  be  approximately  constant. 

Allowance  has  therefore  been  made  in  the  model  for  the  theoretical  state-of- 
discharga  to  be  enhanced  by  chin  amount,  whose  magnitud  '.a  governed  by  an  adjust¬ 
able  parawater  specifying  the  percentage  of  anode  caps,  r  degraded  per  minute  of 
discharge  duration.  The  stace-of-disehargo  of  both  electrodes  is  simultaneously 
increased,  so  it  is  assumed  that  whatever  the  mechanism  for  the  capacity  loss,  it 
involves  lithium  and  iron  disulphide  in  the  same  relative  proportions  as  the 
nominal  electrical  discharge  raacticn  on  the  upper  voltage  plateau.  All  simula¬ 
tions  shown  in  Figs  4  to  14  (except  Fig  !!)  are  based  on  the  same  fixed  capacity 
degradation  rate  of  a  fraction  of  a  per  cent  of  lithium  theoretical  capacity  per 
min.  Experimental  evidence  showed  that  this  wan  the  appropriate  rate  for  one 
particular  battery,  and  it  is  encouraging  that  this  figure  ia  broadly  in  agree- 
ment  with  the  range  recently  quoted  for  isothermal  ai  gle-cells  stood  on  open- 
circuit  before  discharge.  This  range  of  between  2.5  and  10  eoulombs/min  between 
450  and  550°C  for  amall  calls  containing  different  salt  eutectics  equates  to 
0,06  to  0.23  or  0.2  to  0.82  of  lithium  capacity  per  min,  depending  on  whether  the 
whole  of  the  lithium  in  the  anode,  or  only  that  fraction  of  it  associated  with 
the  relevant  constant-potential  phase  transition!  ia  considered, 

No  thermal  effects  of  this  capacity  loss  were  taken  into  account  in  the 
simulations  shown  in  theae  Figures;  this  possibility  is,  however,  considered  in 
section  7.2. 

6 . 2  Allowance  for  time-dependent  polarisation 

The  impetus  for  this  refinement  arose  from  RAE  tBBts  on  Battery  I,  since 
it  was  seen  in  section  5  that  the  Mark  XI  model  simulated  this  battery  poorly 
(Fig  I).  In  this  teat,  the  following  conditions  were  arranged  to  be  met! 

(i)  the  battery  was  fully  discharged  before  cooling  ended  ita  life; 

tii)  the  current  density  was  Bet  so  that  the  resistive  heating  rate  in 

the  cells  balanced  the  natural  cooling  rate  according  to  Newton’s 
law,  us  predicted  by  the  thermal  model,  The  test  thus  approximated 
to  an  isothermal  battery  discharge,  and  this  was  confirmed  by  a 
thermocouple  placed  within  the  insulating  wraps  of  the  battery  near 
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to  the  perimeter  of  the  cell  stack.  One  variable  affecting  perform¬ 
ance  was  thereby  minimised; 

(iii)  the  current  was  interrupted  as  shown  in  Fig  I  so  that  pulse  shapes 

could  ba  analysed. 

The  observed  pulses  have  the  following  features  in  common) 

(i)  when  the  heavier  current  is  applied  over  each  50  a  in  every  60  s  cycle, 
and  then  switched  back  to  the  lighter  load  for  a  further  10  a,  the 
voltage  soon  recovers  to  a  level  only  slightly  below  that  before 
starting  the  pulse  sequence,  especially  after  allowance  is  made  for 
the  increase  in  state-of-discharga  during  the  pulae  cycle. 

(ii)  the  relative  contributions  of  the  instantaneous  versus  the  time- 
dependant  voltage  changes  ars  different  for  tha  recovery  pert  of  each 
60  s  cycle  from  those  for  the  initial  voltage  dip.  Thus,  the  instan¬ 
taneous  recovery  is  larger  than  the  instantaneous  voltage  drop,  whilst 
the  range  of  the  time-dependent  recovery  is  smaller  than  the  time- 
dependent  drop; 

(iii)  the  internal  reaiitance/polarlsation  scarcely  increases  until  the 
battery  is  deeply  discharged,  when  the  on-load  voltage  passes  through 
a  'knee*  before  collapsing! 

(iv)  the  shape  of  both  the  polarisation  and  recovery  time-dependent  pulse 
curves  may  be  approximated  by  exponential  decay  expressions! 

(v)  tests  on  the  same  battary  discharged  under  other  conditions  indicated 
chat  the  magnitude  of  the  time-dependent  polarisation/recovery  was 
approximately  proportional  to  the  applied  current,  as  is  the  instan¬ 
taneous  IR  drop  or  rise. 

The  aim  of  this  version  of  the  modal  has  been  to  search  for  a  mathematical 
expression  having  the  appropriate  form  to  eater  for  continuous  and  pulsed  dis¬ 
charges,  with  continuous  discharges  being  treated  mathematically  as  a  succession 
of  pulses  in  all  of  which  the  current  values  are  identical.  The  expression  should 
incorporate  a  time-dependent  polarisation/recovery  effect  as  well  as  merely  an 
instantaneous  drop  or  rise.  Restrictions  on  the  form  of  mathematical  expression 
which  can  be  used  automatically  result  from  consideration  of  the  following  scen¬ 
arios.  Two  model  simulations  are  run,  whose  only  difference  is  that  in  the  first, 
a  single  long  time  period  is  chosen,  throughout  which  a  single  current  value  (or 
resistor  value)  operates)  in  the  second  simulation,  this  same  time  period  is  sub¬ 
divided  into  two  or  more  segments  in  tha  model  input  data,  in  all  of  which  the 
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Same  current  or  resistor  value  operates  as  in  simulation  I.  E xporience  dictates 
that  identical  results  should  bu  obtained  in  both  cases,  and  that  in  ths  sscond 
case,  the  model  should  not  produce  voltage  discontinuities  at  the  boundaries  bet¬ 
ween  successive  time  periods. 

One  expression  which  obeys  this  restriction  approximately,  though  not 
exactly,  and  yat  gives  an  improved  simulation  over  the  Mark  II  model  1st 

V  -  V*  -  IRf I  +  a( 1  -  e“bt)]  -  I'Rce"dt  ,  (2) 

where  t  ie  the  time  from  the  beginning  of  the  last  load  change, 

R  is  the  internal  resistance  at  a  particular  instant, 

I  ia  the  current  operative  in  the  time  period  in  question, 

I'  is  the  current  operative  in  the  previous  time  period 

and  a,b,c,  and  d  are  constants  to  be  fitted  with  numerical  values. 

To  obey  the  above  mathematical  restriction  it  may  be  Shown  that  a  ■  c  and 
b  *  d  ,  Obviously  I'  is  undafinsd  in  the  first  time  period,  so  in  the  Mark  III 
model  the  second  exponential  term  has  bsen  omitted  for  this  period,  In  physical 
terms  the  expression  provides  an  initantanaoua  drop  as  in  the  Mark  II  moitnl  at 
the  beginning  of  the  first  time  period,  and  the  voltage  then  undergoes  a  time- 
dependent  exponential  decay  at  a  rate  determined  by  the  value  of  b  ,  to  an 
eventual  limit  of 


V  -  V'  -  IR  (t  ■  0)  +  V  -  V'  -  IRC  1  +  a)  (as  t  ttnds  to 

infinity) . 

It  also  reaumss  with  this  value  at  the  atart  of  the  second  time  period  if  the 
current  is  unchanged  and  provided  a  "  c  and  b  ■  d  ,  Provided  the  first  time 
period  is  not  so  short  that  the  voltage  is  still  far  from  its  eventual  limit  by 
the  end  of  it,  the  voltage  will  heve  reached  approximately  this  value  ensuring 
no  significant  discontinuity  at  the  boundary  of  the  first  and  second  time 
periods . 

Whan  ths  current  level  in  Che  second  timo  period  differs  from  that  in  the 
first,  ths  expression  also  provides  for  a  recovery  curve  or  further  polarisation 
curve  which  in  broad  terms  mimicks  ths  actual  behaviour  of  batteries,  This  is 
seen  in  Fig  4,  in  which  curve  1  shows  the  result  of  stepping  up  the  current  to 
increasingly  higher  levels  in  each  successive  60  second  period,  beginning  5  t 
after  activation.  Curve  2  (the  solid  trace)  shows  the  effect  of  alternating  the 
fixed  current  between  two  different  levels  in  these  same  time  periods,  and 
curve  3  (broken  trace)  shows  thst  whan  the  simulation  of  curve  2  is  repeated  only 
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using  die  higher  current  for  all  time  periods  then  a  smooth  prediction  without 

noticeable  discontinuities  is  produced • 

With  the  conetanta  a  equal  to  c  and  b  equal  to  d  ,  equation  (2) 
unfortunately  gives  the  same  voltage  span  during  recovery  as  during  polarisation, 
which  was  not  the  case  for  Battery  I  at  notsd  above,  In  concept,  equation  (2) 
therefore  treats  subsequent  tint  periods  at  a  balance  between  the  polariaing 
affect  of  the  present  currant,  opposed  by  a  relaxation  or  further  polerieation 
from  the  previous  currant  level,  depending  on  whether  it  was  greater  than  or  less 
than  the  present  levsl.  This  notion,  although  crude,  nevertheless  increases  the 
realists  of  this  Mark  III  vtraion  ovsr  tha  Mark  II  varsion  for  nearly  all  the  dis¬ 
charge  data  available  to  RAE,  at  will  bs  shown  balow, 

6 . 3  Improvements  to  model  procedures,  efficiency,  and  output 

The  Mark  II  model  (eection  5)  automatically  produced  e  temperature  distri¬ 
bution  equal  in  site  to  the  chosen  grid  sits  for  each  time  the  reeulta  were 
sampled.  The  estimated  on-load  voltage,  and  eight  other  lingle-valued  functions 
specifying  heat  lost  rates  and  fluxss  were  alio  eamplad  at  each  time.  For 
designs  specif iod  on  large  grids,  and  sampled  many  times,  the  output  is  volumin¬ 
ous.  A  small  selection  programme  is  therefore  useful  to  pick  out  particular 
f  sturts  for  subsequent  plotting. 

In  order  for  tha  Mark  III  model  to  reproduce  curved  pulse  shapes  it  is 
necessary  hoth  to  calculate,  and  to  dieplay  the  reeulta  at  frequent  intervals 
around  the  time  of  each  pulse.  The  provious  method  (section  3,  item  (xvii))  of 
interrupting  the  expanding  tima  intervals  between  calculations  in  order  to 
sample  two  points  either  side  of  e  load  change  would  now  be  inadequate.  Instead, 
the  programme  has  now  been  modified  to  operate  as  follows.  It  remains  as  des¬ 
cribed  in  eection  3  until  the  first  load  change  specified  by  the  user  is 
approached,  If  the  next  calculation  time  would  have  exceeded  the  time  of  the 
load  change,  the  interval  is  reset  so  as  to  make  the  calculation  time  coincide 
with  this  time.  The  following  interval  is  then  reset  to  an  arbitrary  0.3  a , 
Subsequently,  the  calculation  intervals  start  expanding  again  at  the  same  17,  rate 
each  time,  but  from  tha  reset  time  of  0.5  s,  rather  than  the  point  which  they 
had  expanded  to  prior  to  tha  onset  of  the  load  change.  For  pulsss  of  less  than 
1  •  in  duration,  the  0.5  s  reset  level  would  naad  reducing  further.  Examples  of 
the  type  of  model  output  produced  were  shown  in  Fig  4,  end  are  also  seen  in  sub¬ 
sequent  Figures, 

This  modification  has  two  effects:  it  enables  curvature  to  be  reproduced  on 
all  but  the  briefest  of  pulses,  but  to  do  thie,  the  programme  needs  to  perform 
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more  iterations  to  achieve  this  greater  realism.  As  the  overall  computation  time 

increases  approximately  in  proportion  to  tha  number  of  iteration!  performed  (and 
doat  not  aignif icantly  increase  due  to  improvements  in  tha  physic*  during  any  one 
iteration),  eimulaeiona  involving  larg*  numbers  of  pulies  now  take  longer  to  run 
than  previously.  Evan  ao,  all  simulations  shown  in  tha  Figures  were  accomplished 
within  a  period  of  20  to  600  s  using  the  Cray  2  supercomputer  at  RAEj  a  typical 
tim*  would  ba  100  to  200  s  and  only  tha  long-duration  multi-pulse  discharge  of 

battery  3  in  Figs  3  and  9  took  600  s.  Consideration  is  now  being  given  to  making 

1  3 

thia  Fortran  programme  available  for  use  on  advanced  desktop  computers  . 

With  the  Mark  III  modal  now  calculating  and  sampling  the  results  at  many 
more  time*,  thara  was  a  danger  of  breaching  data  storage  limits  if  temperature 
diatributiona  ware  also  bsing  displayed  at  all  these  extra  times.  Also,  whilst 
there  was  a  need  to  calculate  voltages  vary  frequently  around  pulaaa,  there  waa 
no  need  to  dieplay  temperatures  io  frequently.  Hence  the  Mark  III  model  provides 
a  user-chosan  option  to  either  produce  output  in  the  form  of  the  Mark  II  model 
(with  temperatures)  or  else  to  msraly  produce  a  file  of  times  and  voltages  suit¬ 
able  for  direct  plotting,  and  which  can  contain  thousands  of  data  points, 

6,4  Effect  of  discharge  state  on  internal  resistance 

In  common  with  other  types  of  battery,  thermal  batteries  tend  to  exhibit 
a  alight  increase  in  internal  resistance  with  atata-of-discharga ,  aa  shown  in 
Fig  1.  Argonne  and  Sandie  identified  tha  cathode  aa  contributing  moat  to  this 
increase9,  and  other*6’1*  have  shown  that  tha  bulk  resistivity  of  its  first 
discharge  product,  LiijF*^,  is  about  20  times  that  of  iron  disulphide.  However, 
as  explained  for  ona  particular  diacharg*  (Fig  31b  of  Ref  135,  the  battery  inter¬ 
nal  resistance  may  only  thow  about  a  two-fold  increase  at  its  peak  due  to  the 
inhomogeneity  of  discharge  producing  different  phases  simultaneously,  and  the 
fact  that  tha  subsequent  discharge  product  is  much  more  conductive,  Also,  if  the 
electrochemical  balance  between  the  electrodes  is  such  as  to  totally  deplete  the 
lithium,  for  example  in  immobilised  liquid-lithium  anodes,  as  the  end-of-life 
approaches,  this  could  be  another  factor  increasing  internal  resistance. 

ft  waa  shown,  using  computer  fitting  routines,  that  ths  experimental  pulse 

shapes  from  Battery  I  could  best  be  fitted  to  an  expression  involving  the  sum  of 

1 6 

both  exponentiel  and  linear  tarms  .  This  had  previously  been  employed  by 
a 

Shepherd  .  However,  this  form  of  expression  does  net  obey  the  mathematical  res¬ 
triction  noted  in  section  6.2,  since  the  linear  term  causes  the  voltage  to  drop 
(or  rise),  without  limit.  An  alternetive,  which  uees  the  form  of  aquation  (2), 
is  therefors  to  incorporata  the  required  variation  within  the  internal  resistance 
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function  R.  Several  expressions  defining  the  variation  of  the  uempercture- 

compeneated  internal  resistance  with  etate-of-diacharge  were  therefore  evaluated 
against  the  body  of  experimental  discharge  data.  In  one  case,  quite  a  good  fit 
to  the  discharge  of  Battery  1  could  be  obtained  by  combining  equation  (2)  with 
equation  (3) 

R  ■  R'(l  +  1.3Q2)  ,  (3) 

where  R'  is  the  temparature-oorrectud  internal  resistance  at  the  start  of  dis¬ 
charge,  and  Q  ie  the  atata-of-diecharge  of  lithium  in  the  anode  (Q  »  0  at  the 
■tart  of  diechargi)  Q  ■  I  for  total  lithium  exhaustion).  This  ia  aten  in  Rig  3. 
Whan  equation  (3)  was  eubsequently  applied  to  e  number  of  discharges  of  other 
batteries  (Rig  6),  e  consistent  discrepancy  became  apparent,  though  comparison 
of  this  Figure  with  Fig  2  show*  a  nat  improvement  using  the  Mark  III  model  over 
the  Mark  11  veraion  for  most  discharges.  This  was  that  tha  predicted  on-load 
voltage  waa  too  low  at  low  stataa-of-diecharga,  but  then  crossed  over  the  experi¬ 
mental  curve  in  mid-life,  and  finally  bacame  too  optimistic  near  the  end  of  life. 

The  specific  resistivity  allocated  for  the  electrolyte  pellets  (section  3, 
item  (xiii) )  wee  then  updated  with  the  figure  of  <v1.0n  cm  more  recently  quoted 
by  Sandie  I,aboratoriai!  2  for  LlCl/KCl/MgO  pallet*.  Simultaneously,  other 
equations  for  the  instantaneous  internal  resistance  were  evaluated.  One  expres¬ 
sion  which  better  fitted  the  discharges  of  aeveral  batteriaa  subjected  to  con¬ 
tinuous  loads  waa 

R  ■  R'(l  +  2Q  +  Q2)  .  (A) 

This  unfortunately  geve  a  poorer  aimuletion  for  Battery  1  than  that  ihown  in 
Fig  5,  so  it  appears  thet  this  feature  may  need  optimising  for  each  category  of 
battery.  Ae  a  reasonable  compromise,  Bnd  also  to  enaure  that  tha  battery  voltage 
collapses  when  the  lithium  in  the  anode  ia  exhausted,  the  form  of  equation  used 

g 

by  Shepherd  and  others  waa  than  adopted 


This  automatically  provide*  for  a  final  'knee'  in  the  discharge  curves  seen  in 
Battery  1  and  aoua  other  batteries  a*  the  aotive  materials  become  exha  '.ated,  a 
feature  which  equations  (3)  and  (A)  fail  to  provide.  Figs  7  to  14  show  simula¬ 
tion*  produead  by  the  combination  of  aquation*  (2)  and  (5)  and  with  a  •  c  «  0,7 
and  b  *  d  ■  0.2  ,  for  tha  wide  range  of  battery  end  teat  conditions  noted  above. 
Comparison  of  Fig  8  with  Riga  2  and  6,  and  of  Fig  7  with  Fig  1  shows  that  the 
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combination  of  aquation*  (2)  and  (5)  in  the  Mark  ill  (nodal  afford*  a  furthar 
improvement  for  most  of  batteries  !  to  6.  In  no  case  has  this  combination  sig¬ 
nificantly  worsened  the  simulation,  though  comparison  of  Fig  9  with  Fig  3  shows 
little  difference  for  two  batteries  subjected  to  regular  short  pulses  on  top  of 
a  low  background  currant. 

6.5  Effect  of  manufacturing  tolerances 

Both  electrode  active  material*  can  exhibit  different  discharge  plateaux. 
This  is  particularly  apparent  whan  lithium-silicon/iron  disulphide  batteries  are 
heavily  diachargad,  as  aaan  in  Fig  10.  Th*  ainsulation  curvaa  for  different 
external  circuit  raeiitora  era  seen  to  exhibit  sharp  'corners'  whereas  the  cor¬ 
responding  experimental  curves  are  always  smoother.  In  the  model,  all  cells  are 
reaching  identical  stataa-of-diecharga  at  the  same  time,  which  hae  generated  the 
discontinuities, 

Bernard!  and  Newman^  accounted  for  the  rounding  of  the  discharge  curves  of 
lithium-silicon/iron  disulphide  cells  in  terms  of  thermodynamic  irreversibilities 
aeeociatad  with  ohmic  losses,  migration,  msas-tranefar,  and  ovarpotantials.  For 
multi-cell  batteries,  it  is  also  possible  that  tha  greater  smoothness  of  the 
experimental  dischargee  could  be  at  least  partly  dua  to  th*  known  scatter  in  the 
weights,  and  hence  theoretical  capacities,  of  the  active  material  pellota,  due  to 
manufacturing  tolerances  during  tha  cold-preasing  operation. 

A  new  adjustable  parameter  was  therefor*  provided  specifying  tha  range  of 
this  tolaranc*  band  ns  a  fraction  of  th*  mean  weight.  Than,  whan  the  electrode 
theoretical  capacities  had  bean  defined  in  tha  programme,  a  population  of  differ¬ 
ent  capacitiaa  within  this  range  was  also  defined.  The  definition  of  battery 
thermodynamic  emf  was  than  based  on  a  calculation  of  the  equilibrium  potential  of 
each  individual  cell  in  a  aeries  string,  with  its  unique  theoretical  capacity. 
Initially,  a  non-random  distribution  was  sat  up  in  which  the  allowable  range  was 
equally  divided  between  rhe  required  number  of  cells,  with  each  call  assigned  a 
weight  at  a  different  point  on  this  regular  distribution,  Subsequently,  a 
random  number  generator  routine  was  called  by  the  programme  to  provide  a  com¬ 
pletely  random  distribution  in  this  range. 

The  effect  of  this  latest  modification  on  the  simulation  of  Battery  B  may  be 
■ecn  by  comparing  Fig  II  with  Fig  !0|  a  closer  match  has  been  achieved  by  incor¬ 
poration  of  t  very  real  effect.  The  simulations  in  Fig  II  were  generated  by 
uomewhat  reducing  the  capacity  degradation  rate  from  the  level  uted  in  all  other 
simulations,  and  also  by  assigning  s  maximum  spread  of  pallet  weights  of  51 
either  Bide  of  tha  mean.  (Thie  was  done  for  illustration  only  and  is  not  neces¬ 
sarily  representative  of  actual  prartice.)  The  simulations  shown  in  Figs  12  to  14 


TM  MS  1163 


18 

also  used  this  same  weight  distribution,  hut  reverted  to  the  capacity  degradation 
rata  previously  uaad.  Fig  12  show*  that  for  Battery  JO  in  particular,  the  MarktII 
modal  incorporating  thia  effect  broadly  reproduce!  the  inflexions  aean  in  the 
experimental  diecharge  curve,  a*  far  ae  it  is  possible  to  measure  them  by  ruler 
from  the  report  (Fig  31a  of  Ref  15). 

7  TEMPERATURE-LIMITED  BATTERY  DISCHARGES 

There  is  no  doubt  that  the  voltage  fall-off  observed  in  some  discharges  may 
be  linked  to  the  freesing  of  molton  salts  rather  than  to  active  materiel  exhaus¬ 
tion.  Table  I  shows  which  dischargee  ware  predicted  to  be  temperature-limited 
within  the  tlmeacale  of  tha  varioua  figurai.  In  improving  the  voltage-time  model 
further  for  such  discharges,  it  is  aaaantlal  to  accurately  predict  the  elapsed 
time  before  molten  salt  in  tha  battery  re-solidif ies,  The  prediction  of  this  time 
may  be  affected  by  uncertainties  in  two  independent  factors  -  the  freecing  tempera¬ 
ture  and  the  elope  of  the  cooling  curve  -  both  of  which  are  discussed  below. 

7 . 1  Salt  freesing  temperature  in  discharged  batteries 

There  is  good  agreement,  within  a  few  degrees  Celeiue,  about  the  freezing 
point  of  pure  salt  eutectics  such  ae  lithium  chlorida-potasaium  chloride,  How- 
evsr  it  is  not  completely  clear  at  what  temperature  molten  aalte  freeze  in  dis¬ 
charged  thermal  batteries.  Experimental  evidence  using  thermocouples  has  in  some 
cases  shown  a  freezing-point  arreat  in  the  cooling  curves  of  cell*  very  close  to 
the  theoretical  point  of  352°C  (see  for  example  Fig  12  of  Ref  18).  A  freezing 
exotherm  near  to  this  theoretical  temperature  wee  also  seen  in  thermocouple  tracts 
for  the  cooling  of  electrically-inactive  'dummy'  bstteriee  containing  LiCl/KCl/MgO 
and  iron/pocaasium  perchlorate  pyrotechnic  a*  the  only  stack  components^.  How¬ 
ever,  thermocouple  measurements  from  tha  same  study  but  involving  elactrienlly- 
active  batteries  suggested  that  freezing  could  occur  over  a  range  of  temperatures 
below  the  theoretical  point.  Contrarily,  Searcy  and  Armijo  have  pointed  out  how 
the  fraaeing  temperature  of  LiCl/KCl  could  be  raised  to  390  to  420°C  in  a  polar¬ 
ised  tharmal  battery,  depending  on  tha  extent  to  which  concentration  gradients 

£. 

build  up  in  it  by  the  passage  of  different  current  level*,  According  to  Newman  , 
this  same  affect  may  lead  to  tha  predicted  precipitation  of  the  salt  components 
within  the  electrode*  of  flooded  celle. 

Certainly,  solidification  can  only  occur  at  a  higher  temperature  if  the 
relative  pioportiona  of  LiCl  to  KC1  change  without  any  contamination  of  the  liquid 
by  ions  not  already  present  in  it.  However  the  possibility  also  exists  that  ions 
such  as  sulphide  or  oxide,  which  are  known  to  have  some  solubility  in  molten 
LiCl/KCl,  may  dissolve  in  it  in  heavily  discharged  batteries  end  have  s  depressant 
effect  on  solidification  temperature. 
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Basad  on  the  apparent  freezing  point  dapreaaion  seen  in  earlier  battery 

4 

tasta  using  internal  thermocouples  ,  all  simulation  curve*  shown  in  the  Figures 
have  bean  generated  uting  302-322°C  as  the  notional  freesing  rang*  of  molten 
ealt.  However,  since  Table  I  shows  that  not  all  of  the  discharges  were  limited 
by  temperature,  it  was  expected  (and  confirmed  by  further  model  simulations  using 
the  theoretical  freesing  point  of  352°c)  that  not  all  simulation*  would  be  affec¬ 
ted  by  a  change  in  assigned  frsecing  point.  Referring  to  Figs  7  and  S,  the  pre¬ 
dicted  discharges  for  Battsriss  I  and  4  remained  unaffected  by  an  increase  in 
freesing  point  from  302-322  to  332°C|  that  for  Bathary  2  improved  even  more  after 
400  « |  that  for  Battery  6  vaa  a  definite  improvement;  whilst  simulations  for 
Batteries  3  and  3  ware  somewhat  woraened. 

The  two  dlachargea  of  Battery  I!  seen  in  Fig  13  illuatrare  how  the  voltage- 
time  model  can  be  effected  by  eaaigned  freesing  temperature.  When  the  simula¬ 
tions  generated  in  Fig  13  using  a  fretting  range  of  322-302°C  are  compered  with 
those  of  Fig  14  generated  using  the  theoretical  freesing  point  of  332°C,  it  is 
assn  that  the  predicted  discharge  at  the  maximum  ambient  temperature  remains 
unaffected  by  the  change,  indicating  thee  the  cell  temperature  in  the  model  hat 
not  reached  either  framing  noint.  In  contract,  an  inerseee  in  easigned  fretting 
point  haa  improved  the  aitnulation  of  the  discharge  curve  at  the  minimum  ambient 
temperature ,  indicating  chat  duration  wee  predicted  to  be  limited  by  cell 
h'empereture  in  thi*  cue,  It  is  alsu  noteworthy  how  well  the  model  is  able  to 
predict  the  effect  of  ambient  temperature  on  battery  voltage  earlier  in  the  dis¬ 
charge,  though  after  voltage  'spiking'  haa  ended. 

On  balance ,  and  in  thu  absence  of  any  allowance  for  the  other  heat  aouraea 
mentioned  in  the  next  section,  an  assigned  fvessi.ng  point  at  or  below  the  theo¬ 
retical  freesing  point  of  (but  not  above  it)  eeeme  beet  for  obtaining  tha 

beet  prediction  of  battery  duration  in  temperature-limited  discharges.  However, 
further  ntudy  of  thia  aspect  would  be  desirable. 

7 . ?  Other  heat  source! 

Cooling  rates  of  thermal  batteries  are  deelgned  to  be  email,  *o  their  cool¬ 
ing  curves  will  automatically  intersect  any  cunetant-temparature  line  at  a  glanc¬ 
ing  angle,  to  magnifying  tha  uncertainty  of  when  the  voltage  is  predicted  to 
oollapea  if  the  electrolyte  freasing  temperature  is  not  precisely  known.  Another 
source  of  uncertainty  arises  if  the  elope  of  the  cooling  curve  should  bs  affected 
by  extra  heat  aourcaa.  In  deriving  tha  aimvlationa  illustrated,  the  only  heat 
aourcei  takati  into  account  were  from  the  pyrotechnic  burning,  heating  due  to  cur¬ 
rant  flow  inside  the  battery,  and  entropic  cooling. 
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Rigorous  theoretical  treatments  of  heat  sources  within  batteries  have  been 

given  '  which  includa  a  few  effects  not  preiently  modaliad.  However,  it  wat 
pointed  out  that  some  of  theta  could  be  negligible  for  tome  type*  of  battery.  In 
any  cate,  theta  treatment*  deal  only  with  affteta  associated  with  the  theoretical 
operation  of  the  battery  rather  than  with  othar  unintended  eventualities  auch  as 
leakage  currents  or  cathode  pyrolytit.  These  letter  may  ponibly  dwarf  aome  of 
the  theoretical  af facta  during  thermal  battery  discharge*. 

tiling  curve-fitting  procedures  it  was  estimated  that  up  to  38  cal  of  heat 
could  be  produced  per  gm  of  cell  in  a  csicium/calcium  chromate  battery  ,  rep¬ 
resenting  ebout  30X  of  the  hast  output  from  tha  pyrotechnic.  However,  for  a 
magnasium/iron  disulphide  battery,  no  extra  heat  had  to  be  invoked  to  fit  experi¬ 
mental  data.  Examples  have  baan  given4  of  lithium/ire.i  disulphide  batteries  whoso 
temperature*  could  also  be  modelled  satisfactorily  without  recourse  to  extra  heat 
generation. 

Despltu  this,  the  rare  occurrences  of  catastrophic  failures  of  thermal 
bacteria*  are  a  reminder  Chat  tha  energy  of  the  slactroactive  materials  is  com- 
parabls  with  that  of  the  pyrotechnic.  A  consideration  of  their  thantiodynamica, 
their  elactrical  efficiency,  end  the  rates  of  side  reactione  tharefor*  helps  to 
bracket  tha  magnitude  of  any  temperature-time  effect*.  The  Appendix  indicates 
one  way  in  which  tha  scale  of  hast  production  from  chemical  aide-reactions  majt  be 
•  stimatad.  This  fit  at  consider*  tha  total  heat  of  chemical  reaction  if  ell  the 
lithium  metal  were  depleted)  it  then  eitimatt*,  for  the  degradation  rate  of 
electrical  capacity  mentioned  in  section  6.1,  how  this  scaled-down  heat  produc¬ 
tion  rate  might  affect  the  cooling  curve*  of  tha  different  batteries  in  Table  I. 
This  analysis  does  not  identify  tha  mechanism  of  haat  production,  but  linking  it: 

to  the  linear  degradation  rat*  discussed  in  section  6,1  would  be  consistent  with 

1  2 

either  the  continuoua,  slow  direct  chemical  reaction  recently  postulated 
between  dissolved  lithium  snd  iron  disulphide,  or  with  slight  but  continuous 
shorting.  It  would  not  be  so  consistent  with  cathode  pyrolysis  which  might  be 
expected  to  peak  early  in  the  battery's  active  life, 

The  effect*  of  this  heat  production  rale  were  investigated  for  all  of  the 
batteries  in  Table  I.  For  the  small  Batteries  3,  3  and  6,  it  made  an  insignific¬ 
ant  difference  to  the  predicted  cooling  curve,  Tha  cooling  ratea  of  Batteries  1, 

2  and  4  were  slightly  reduced,  but  those  for  Batteries  8  to  1 1  were  significantly 
reduced  ai  a  result  of  chemical  heat  generation,  The  comparisons  for  three 
butteries  are  shown  in  Fig  15.  The  larger  affects  for  Batteries  S  to  1 1  could  bo 
due  toi 


TM  MS  1163 


2 1 


(i)  tha  low  #urface/voluu»  ratio  for  large  r.Blla,  reducing  heat  loaa; 

(ii)  the  fact  that  vary  good  inaulation  waa  uaad,  becauae  of  the  long 
duration  requirements ,  which  further  reduced  heat  least 

(iii)  the  fact  that  these  lithium  alloy  batteries  contained  a  higher 
weight  percentage  of  lithium  per  cell,  because  of  the  desirability 
of  discharging  the  battery  on  juac  one  of  the  voltage  plateaux  of 
the  alloy. 

Of  course,  these  calculations  era  for  illustration  only,  and  they  might 
need  to  be  modified  if  any  of  the  assumptions  made  proved  to  be  unwarranted. 

For  example,  not  all  of  tha  lithium  might  be  available  to  engage  in  axothermit 
reaction*,  or  conversely  other  reducing  agents  like  silicon,  aluminium  or 
powdered  Iron  might  also  contribute  to  extra  heat  production. 

It  is  seen  from  Fig  15  that  the  predicted  time  elapsed  before  aome  bet- 
teries  (especially  Battery  II)  reached  a  fixed  temperature  would  be  longer  where 
there  was  an  extra  heat  aourca  than  without  it.  Alio,  when  this  effect  it  linked 
to  freeeing  temperature  (section  7.1)  in  order  to  predict  the  point  at  which 
voltage  fall-off  is  expected,  the  freeeing  temperature  in  the  model  may  need  to 
be  set  higher  when  thia  haat  aourca  is  incorporated  than  without  it, 

fl  REMAINING  DISCREPANCIES 

Figa  6  to  14  show  that  the  latest  improvements  made  in  the  Mark  III  model 
have  provided  a  net  increase  in  the  realism  of  the  eimulations  for  nearly  all 
the  battariee  and  teat  conditions  which  RAE  possess  information  on.  Even  so, 
inspection  of  the  Figures  highlights  the  following  features  of  ths  simulations 
requiring  further  improvement! 

(i)  the  simulated  durations  may  still  be  Coo  optimistic  for  batteries 
which  experience  very  high  temperature#  over  a  coneidarable  timet 

(ii)  for  batteries  subjected  to  occasional  pulses  with  open-circuit 
periods  in  between,  the  simulated  voltage  between  each  pulse  does 
not  fall  away  as  it  did  in  practice  (Batteries  3  and  7,  Fig  9), 

Also,  the  simulated  voltage  drops  during  the  pulses  do  increase  as 
the  discharge  proceeds,  as  observed,  but  the  increases  are  too 
sudden  and  too  large.  Further  study  of  the  temperature-dependence 
of  the  electrical  resiatance  of  solidified  electrolyte  pellets  might 
help  to  improve  thisi 

(iii)  whilet  the  treatment  of  time-dependent  polarisation  following  appli¬ 
cation  of  a  heavier  load  is  adequate,  the  simulated  magnitude  of  tha 
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recovery  curve  efter  the  load  ha*  been  removed  1b  too  large 
(Battery  I). 

9  CONCLUSIONS 

(i)  A  further  stage  In  the  refinement  of  eh*  RAE  voltage-time  model  of 
thermal  batteries  ha*  resulted  in  greater  realism  of  the  resulting  simulations, 
over  those  produaad  by  the  previous  version. 

(ii)  Application  of  the  naw  model  to  a  limited  number  but  a  wide  ranRe  of 
battery  types  and  tsse  condition*  hat  shown  that  both  the  overall  voltage  level 
and  the  internal  resistance  of  the  majority  of  batteries  are  estimated  quite 
well,  especially  early  in  the  discharge  life,  but  diverge  from  the  observed  data 
latar  in  life. 

(iii)  iulse  shape*  may  now  be  simulated  batter,  though  further  improvement 
is  neceeeary. 

(iv)  If  the  whole  of  the  lithium  in  the  anode  (but  no  other  reducing 
agent)  is  assumed  to  be  available  to  engage  in  chemical  aide  reactions,  and  the 
lithium  is  being  depleted  by  this  means  at  a  rate  of  a  fraction  of  s  per  cent  per 
minute  of  active  life,  it  ie  predicted  that  tha  haat  rtlane*  would  searcaly 
reduce  the  cooling  rates  of  many  small  to  medium-sised  batteries,  However,  it 
we*  predicted  to  vary  considerably  reduce  tha  cooling  rate  of  eoma  large,  long- 
duration,  well-insulated  battariaa  containing  lithium-alloy  anodee.  This  might 
delay  the  point  at  which  their  voltages  fell  away  dua  to  freaaing  of  molten  salts. 

(v)  Whilet  o  reasonable  simulation  of  performance  over  the  whole  range 
of  battery  type*  may  be  obtained  by  just  ons  set  of  equations  end  parameter 
values,  there  ere  indication*  that  further  optimisation  of  the  simulation  for 
particular  categories  (such  as  the  products  of  one  particular  manufacturer,  or 
one  particular  combination  of  materials)  may  be  possible  by  resetting  parameter 
values  individually  for  that  category.  Better  simulations  may  then  result  when 
the  model  ia  applied  to  new  battery  applications  within  that  category.  A  uni¬ 
versal  framework  has  been  developed  within  which  to  do  thii  -  it  would  merely 
involve  readjustment  of  numerical  parameters  based  on  discharge  curves  for  o  few 
typical  batteries  in  each  category. 

(vl)  The  model  is  not  Intended  to  eimulete  the  performance  of  those  bat¬ 
teries  which,  through  some  manufacturing  defect  or  other  meana,  generate  excess¬ 
ive  heat  or  develop  serious  short-circuits.  However,  a  steady  continuous  inter¬ 
nal  short  may  be  simulated  to  predict  its  effect  on  both  battery  temperature  and 
voltage. 
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Appandix 

ESTIMATED  MAGNITUDE  OF  HEAT  GENERATION  FROM  SIDE  REACTIONS 
(see  section  7.2V  " 

As  s  first:  approximation,  it  may  bo  assumed  that  exotherms  are  generated 
by  direct  chemical  reaction  of  lithium  with  either  iron  disulphide  or  with  other 
materials  such  us  oxides  in  the  thermal  insulation.  Interaction  with  iron  disul- 
phide  could  either  be  by  direct  physical  contact,  or  via  pyrolysis  generating 
sulphur  vapour  and  migration  of  thia  to  the  anode.  It  could  also  occur  by  dis¬ 
solution  of  tha  materials  in  elactrolyts  molten  salts,  and  subsequent  diffusion 

1 2 

of  these  into  Che  separator,  as  recently  suggested  .  Whatever  the  mechanism  it 
is  assumed  that  tha  heat  release  is  the  same  as  that  calculated  for  the  nominal 
discharge  reaction,  though  brought  about  entirely  by  chemical  means,  This  cor¬ 
responds  to  the  case  where  no  useful  electrical  work  is  done,  and  where  all  the 
energy  ie  wasted  as  internal  heat. 

The  following  treatment  will  assume  an  immobilised  liquid-lithium  anode 
in  an  electrically-balanced  cell,  according  Co  the  nominal  chemical  or  electro¬ 
chemical  rsaction 

3Li  +  ?KeS2  +  Li^e^  . 

If  the  heat  release  from  thie  may  be  approximated  by  the  change  in  enthalpy  at 
any  temperature  T  ,  then  literature  values  for  tha  cell  emf  and  its  Lamporature 
coefficient  may  be  used  in  the  Dlbba-Helraholr.  equation  to  calculate  the  heat 
generation  expected  from  complete  conversion  of  the  amounts  shown  by  the  above 
equation.  This  yields  a  value  of  130  kcal,  or  *■6250  cal/gm  of  lithium  metal. 

For  a  cell  containing  about  2.25  wtX  of  lithium,  the  output  would  be 
~U0cal/g  of  cell  if  all  tha  lithium  metal  became  exhausted  by  unproductive 
reaction.  Assuming  a  weight-average  cell  specific  heat  of  about  0.2  cal/g/°C, 
and  that  all  the  active  ingredients  could  be  brought  together  instantly,  the 
'battery'  would  then  approach  a  bomb  calorimeter  in  behaviour,  with  an  instant 
temperature  rise  of  -700°C  over  and  above  that  produced  by  tha  pyrotechnic. 

When  a  simulation  was  run  for  one  of  the  batteries  using  the  short-circuit  option 
mentioned  in  section  5,  and  resistances  of  I  ME  and  1  mfl  were  assigned  for  the 
external  circuit  and  a  continuous  short  down  the  bettery  central  axis  respectively, 
the  temperature  was  predicted  to  rise  by  about  660°C  after  I  min.  This  is  in 
line  with  the  temperature  rise  calculated  above. 

With  regard  to  the  observed  slow  degradation  of  capacity  mentioned  in 
section  6,1,  this  same  heat  output  may  be  scaled  down  in  the  model  (a)  to  that 
fraction  of  the  battery's  capacity  which  is  degrading,  and  (b)  spread  out  over  the 
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Table  I 

DETAILS  OF  BATTERIES  AND  TEST  CONDITIONS 


Battery 

number 

Manufacturer 

Figur® 

refs 

Anode 

type 

Outside 

diameter 

(mm) 

Load 

type 

Predicted  factor 
limiting  life 

HH 

MSA  (Britain) 

1,5,7.15 

Li 

89 

a* 

Capacity 

B1 

SAFT  (UK) 

2,6,8,15 

LiAl 

62 

b 

Marginal 

Kp 

MSA  (Britain) 

2.6,8 

3,9 

Li 

42 

C 

tt* 

Temperature 

Temperature 

■ 

SAFT  (UK) 

2,6,8 

LiAl 

32 

b 

Marginal 

5 

SAFT  (UK) 

2,6,8 

LiAl 

32 

b 

Temperature 

6 

MSA  (Britain) 

2,6,8 

Li 

18 

c** 

Temperature 

7 

MSA  (Britain) 

2,9 

Li 

50 

a 

Temperature 

8 

Sandia  Labs 

10,11 

LiSi 

73  . 

C 

Capacity 

Sandia  Labs 

12 

LiSi 

76 

c 

Temperature 

Sandia  Labs 

12 

LiSi 

51 

c 

Temperature 

H  “ 

Sandia  Labs 

13,14,15 

LiSi 

91 

c 

Capacity  (70°C) 
Temp  ( -35°C) 

Key i  a,  constant  current*  with  intermittent  pulses 
t,  continuous  constant  currant(s) 
c,  continuous  constant  resistive  load(s) 

*,  tested  by  RAE;  the  teat  on  Battery  3  used  b  special  thermal  insulating 
jacket  around  the  battery  to  minimise  heat  loss 

**,  the  figures  Bhow  discharges  at  both  ambient  tempernturo  extremes 

Li,  this  represents  unalloyed  immobilised  liquid  lithium 


Literature  references  to  Sandia  Laboratories  batteries: 

Battery  8:  Kig  8  of  Ref  18 

Battery  9i  Battery  228  in  Fig  5  of  Ref  19 

Battery  10:  Fig  31a  of  Ref  15 

Battery  111  Fig  4  of  Ref  22 
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